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We introduce the anisotropic two-photon Rabi model in which the rotating and counter rotating 
terms enters along with two different coupling constants. Eigenvalues and eigenvectors are studied 
with exact means. We employ a variation of the Braak method based on Bogolubov rotation 
of the underlying su(l, 1) Lie algebra. Accordingly, the spectrum is provided by the analytical 
properties of a suitable meromorphic function. Our formalism applies to the two-modes Rabi model 
as well, sharing the same algebraic structure of the two-photon model. Through the analysis of the 
spectrum, we discover that the model displays close analogies to many-body systems undergoing 
quantum phase transitions. 
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I. INTRODUCTION 

Two-photon Rabi type models serve the study of var¬ 
ious quantum effects in systems of bosonic fields cou¬ 
pled to a set of localized levels. Its simplest instance 
traces back to the theory of micromaser [3, [2] . In such 
a quantum optics set up, three atomic levels are coupled 
each other through two-boson field in a cascade transi¬ 
tion \g) —1 |i) —1 |e), with the energy difference coeg set 
to twice of boson field frequency oj, and with the inter¬ 
mediate state being strongly detuned with Uei and Uig. 
After adiabatically eliminating the intermediate state, 
and neglecting the Stark shift, one arrives at the two- 
photon Rabi model Q. Phenomenologically, such type of 
models describe two-level atom interacting with squeezed 
light [3|. More recently, models of these type emerged in 
different context of quantum technology. Indeed, they 
can be obtained as effective models describing quantum 
dots inserted in QED micro cavity [^; there, the localized 
levels are provided by the excitonic energies of the quan¬ 
tum dot; two photon process in quantum dots were ob¬ 
served through photoluminiscence experiments The 
two-mode Rabi models also emerge in circuit QED in¬ 
volving superconducting qubits in the ultra strong regime 
where non-linear couplings become realistic (7H9| and in 
ion traps 0- Finally, the two-mode boson field can be 
realized by a charged particle in a magnetic field 

In a weak coupling limit between the bosonic degree 
of freedom and the localized levels (that can be modelled 
by a spin degree of freedom), the features of the system 
can be obtained employing the so-called Rotating-Wave 
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Approximation (RWA) (lll| . In most, if not all, of the 
scenarios depicted above, however, the ‘physical working 
point’ is in a strong coupling regime, where the RWA is 
a poor approximation. The limits of the RWA for the 
two-boson Rabi models were investigated numerically in 
0 . 

Here, we deal with the exact energy spectrum and 
eigenstates of the two-photon Rabi model. The exact 
solvability beyond RWA was investigated by means of 
Algebraic Bethe Ansatz in [l^. More recently, the ex¬ 
act eigenenergies and eigenstates were obtained in the 
‘isotropic’ case where the rotating terms and counter¬ 
rotating terms have the same weight [Ii-0. The so¬ 
lution is obtained applying a variation of the procedure 
that Braak recently devised for the single photon Rabi 
model 0 (see also 00). 

In this paper, we introduce a class of anisotropic ver¬ 
sion of the two-mode and two-photon Rabi model, where 
the rotating and counter rotating terms play along with 
two different parameters. As the single-mode case, our 
model enjoys a Z 2 parity symmetry. Relying on it, we 
will give the exact solution of anisotropic two-photon and 
two-mode Rabi model by Bogolubov operator method 
0 0. We will demonstrate that such a model the 
energy levels changes a function of a suitable control pa¬ 
rameter with a similar effects displayed in quantum phase 
transitions of many-body systems. A discontinuous tran¬ 
sition is controlled by the anisotropic parameter. For 
higher values of the spin-boson strength, we find a con¬ 
densation of levels breaking implying a super-radiance 
phenomenon. 

The paper is organized as follows. In Sec. B, the math¬ 
ematical implications arising from the relation between 
the sm(1, 1) Lie algebra and the two-photon Rabi model is 
discussed. In Sec. IB, we present our approach in detail 
for the anisotropic two-photon Rabi model. In Sec. IV, 



2 


we discuss the kind of ’critical behavior’ that we found in 
the energy spectra of the model. In Sec. V, we draw our 
conclusions and future directions. In the Appendix A we 
discuss a specific circuit QED realization of the model 
Hamiltonian. In Appendix B, we sketch the application 
of our approach to the two-mode Rabi model. 


II. THE TWO-PHOTON RABI MODEL: su(l, 1) 
AND Z2 SYMMETRY 

The Hamiltonian reads 

H2-ph = uja^a + Aaz + g{a~^a'^ + 

( 1 ) 

where the uj fix the energy of the bosonic field and g is 
the coupling constant, A being the anisotropic param¬ 
eter characterizing the difference between rotating and 
counter-rotating spin-boson interactions. 

The features of the Hamiltonian above are intimately 
connected with sm(1, I) Lie algbra realized by 

Ko = i(a^a -f i), K+ = AT. = (2) 

obeying to 

[Ko, K±] = ±K±, [K+, K_] = -2Ko, (3) 

with invariant Casimir operator 

C = K+K. + Koil - Ko) = k{1 - k), (4) 

where k labels the representations. K± and Ko are off- 
diagonal and diagonal operators in the Cartan basis, re¬ 
spectively. Relations dS]) can be realized through first or¬ 
der differential operators acting on a suitable functional 
space: 


Kq z —— -t- K , K^ 2^—— -t- 2kz , A_ yy —— . (5) 

dz dz dz 

The isomorphism above can be seen as a ‘Bargman re¬ 
alization’ for sm(I, 1) (see for example [Hi). We remark 
that sit(1,1) is a non-compact algebra and therefore the 
representation space is a non trivial manifold. In the 
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su(1,1) specific case, the representation space is isomor¬ 
phic to an hyperboloid. Different representations corre¬ 
spond to a discrete series in the upper or lower branch of 
the hyperboloid, often denoted as . In the case of the 
two-boson case considered here, k = |;, | for even or odd 
number Fock basis, respectively (see the appendix for the 
representation index for the su(l,I) Rabi model realized 
by two-bosonic-mode). Another possible, but inequiva¬ 
lent, representation is the continuous one, labeled by the 
eigenvalues of K± [l^). As we shall see in the next sec¬ 
tions, the computation of the exact solution of the model 
m involves the diagonalization inside the sm(1, I) alge¬ 
bra. We will choose a representation in which Kq is di¬ 
agonal. Correspondingly, certain constraints will emerge 
in the solution (see (1151) 1. In appendix, we discuss the 
spectrum beyond such constraint. 

The Hamiltonian ([T] ) enjoys a Z 2 symmetry: 


n2-ph = 


cos[—a'a) + % sin(—a' a) 


( 6 ) 


whose eigenvalues belong to a four-dimensional manifold 
spanned by p = ±1 and ±L Such a manifold yields 
four irreducible subspaces of the sm(1, 1) for two photons. 
Incidentally, we note that the model with Z 2 symmetry 
breaking, i.e. with ax term added to the Hamiltonian, 
can also be solved by our method [^ . 

The following observation is important for the proce¬ 
dure to analyze the eigen-system of ©• Despite, the 
(unitary) transformation 

—>■ ia^ , a —>• —ia (7) 


is not a symmetry of the Hamiltonian, it defines an au¬ 
tomorphism of su(l,l). In this way, we shall see that 
© induces a specific redundancy of the eigevectors of 
©, that is ultimately useful for us. The issue is more 
evident by spelling-out the spin basis in 0 


(TT\_( uJo^a + A 5(0^-I-Aal'^) \ 

^ ^ “ V wata - A ) 


( 8 ) 


and applying a rotation of the basis: W 


1 

7^ 


1 

1 


-1 

1 


W\H)W 


1 f 2u}aXa + g{\ + —2A — g{X—l){a^ — a^^) A 

2 —2A -I- p(A — l)(a^ — a^^) 2uja^a — g{l + X){a^ + J 


(9) 


By inspection of the equation above, we see that its 
eigen-system 

»'t(ff)M.-(”;)=E(”;) (10) 


is left invariant by the transformation © by simultane¬ 
ously swapping vi and V 2 - 

Summarizing: the sit(l, I) automorphism ([7]) im¬ 
plies the existence of two sets of eigenstates of 
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O: the original ones {vi,V 2 )'^ and {v 2 ,v[)’^ where 
{v'i,V2)^ = {vi,V2)'^. 

Equivalently, the role of ([71) in the Bogolubov scheme 
we will adopt below is played by the transformation 
z —>■ —z in the sm(1, 1)-Bargman realization (O of the 
spectral problem for ([8]). In this way, the redundancy 
property of the eigenvectors involves (vi{z),V 2 {z)) and 
(u 2 (—z), v[{—z))'^. Indeed, the resulting structure of the 
Hilbert space implies that the domain of analyticity of 
the eigenvectors {vi{z),V 2 {z))'^ can be extended to the 
whole complex plane. 

In the next section, we shall see how such a property 
will be exploited to construct a meromorphic function 
(G(z)) whose analytical structure provides the spectrum 


of the Hamiltonian O- 


III. EXACT ANALYSIS OF THE 
EIGENSYSTEM 

In this section, we will provide the exact analysis of 
spectrum & eigenstates of the anisotropic two-mode Rabi 
model Eq.([T]). We will be formulating an ansatz for the 
eigenvectors expressed as a series expansion that we even¬ 
tually determine through recurrence relations. It turns 
out that such recurrence relations are easier to solve if 
we rotate the spin axes suitably: {H)' = U{H)U\ with 
cos (3 — sin ,8 \ 
sin f3 cos j3 J ’ 


^ wa^a + p + r{a^ + —q + sa^ -I- 


where p = A cos 2/3, q = A sin 2/3, r = 


—q + -I- ta^ ua^a — p — r(a^ -|- 

sin. 2i5 

—-—(1 + X)g, s = (cos^ /3 — Asin^ f3)g, and t = (Acos^ /3 — sin^ l3)g. 


( 11 ) 


r 


Next, we perform a rotation in the sm(1, 1) algebra by 
using a Bogolubov transformation, 

a = ua + va^, = ua’+ va, (12) 

where [a, al] = 1 if = 1. Because of the 

anisotropy, we notice that we cannot get rid of the non- 
Cartan generators of su(l, 1) in all the matrix elements. 


Nevertheless, we can fix the parameters u, u, and /3 in 
such a way that the upper (lower) off-diagonal element 
of {H)' is a generalized lowering (raising) operator in 
sit(l, 1). Such a choice will simplify the solution of the 
recurrence relations below (see (IT^ . (l20l) . (1221) 1. The 
Hamiltonian, then, reads 


(H)' = 


u;77(q!1q! + \) — A cos 2/3 
—Asin2/3 — r(l — A) —(2aAa -I- 1) -b (1 — X)ga^^ 


—A sin2/3 — r(l — A)^(2Q;la -I- 1) -b (1 — X)gcx^ 


w(| -g){a'^ot+ ^) - I - Acos2/3- 


-b a^) 


with 


u = 



V = 



cos 2/3 


1-A 

TTa^’ 


(13) 


where 


V = 


1-(1 + A)V/w2 

1-(1-A)252/w2’ 


(14) 


with the condition 


l5l< 


CJ 

|T+^' 


(15) 


The Eock basis for the new bosonic operators a is 
spanned by \ni)a = (q:1)'"|0)q , m = 0,1,2,..., where 
the vacuum q;|0)q, = 0 is 



We observe that, although the normalization coefficient 
can be fixed by the series expansion of l/-\/u, we work 
with non-normalized Fock basis (the normalization con¬ 
stant does not change the results). 

We are now ready to formulate the ansatz for the eigen¬ 
vectors of (H)': 


(H)' 




As announced in SectHTl the constraint m arises be¬ 
cause of the non trivial topology of the su(l, 1) represen¬ 
tation space. 


(17) 
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(1^1 

02 


J2m=0 


^)c 


Sm=0 ^m\m)c 
where Km and Lm are coefficients to be determined. 


(18) 


Plugging (fTSl) into the spectral problem of the Hamil¬ 
tonian (HU, we obtain, 


fmLm\m)a “ dmKm\'m)a + (1 “ A)g m{m - l)Km\m - 2)a = 0, 


m—0 


m=0 


m—0 


- y] dmLm\m)a + {1 - X)g'^ L 

E 


2)c 


(19) 

( 20 ) 


m—O 


m—O 


m—0 ^ 


(2 \ , 1, w 


2r 


Km\m)a - Km [\m + 2)a + m{m - \)\m - 2)a\ = 0 , 

7? 

' m=0 


where dm= g + r{l — X)(2m + 1)— , fm = ujrj(m +^) — ^+p — E and [a, = maX'^ ^ were exploited. 

u) 2 2 I_ 


By inspection of Eq. (HI, we obtain 

dmKm - (1 - A)(to -I- 2)(to -I- l)gKm +2 


Lm. — 


frr 


( 21 ) 


leading to a closed Ea. (l20l) for Km- 

O-mKm+2 — bmKm “t“ CmKm—2: (22) 


, _ _ 

'-'m. — 


dmjl - X)g ^ 

fm V 

dm (1 — Xf'mijin — l)g‘^ 


(m -I- 2)(m -I- 1), 




fm—2 


I - - 




Cm. — 


(1 - X)gdm -2 _ ^ 

fm-2 V ' 


With a similar logic employed in [l7|, [l9|, 12^, we want 
to extract the eigenvalues of L[ by looking at the ana¬ 
lytical structure of a meromorphic function G{z), con¬ 
structed by imposing that the eigenvectors of the Hamil¬ 
tonian are indeed analytic in whole complex plane. In 
our scheme, G{z) is constructed resorting the property 
of the eigen-system under the application of the sit(l, 1) 
automorphism 0. Namely, the transformation © on 
the new bosonic operators a, of leaves the Hamiltonian’s 
spectrum invariant, with eigenvectors changing as 


^1 

^2 


= G 


P2 




Resorting to the parity symmetry operator Yi 2 -ph , in the 
even Fock space of a, the constant C = ±1, while it is 
C = ±f, for Fock states with odd parity. The transcen¬ 
dental functions in even and odd parity sectors are 

G\^± = (0|[v 52 T V^i], (24) 

= (1| [—'*(^2 T (25) 


where G^ is defined real by multiplying an unimpor¬ 
tant factor —i. Indeed, even and odd functions can be 
considered on equal footing multiplying by the overall 


factor Dm, G' 


A,± 


superscripts ‘o’ and ‘e’, we have. 


G^ . Dropping the even and odd 


G^^-)- — ^ ( ( cos 13Lm 3~ sin[3Km)Dm: 

m—O 

CO 

Gx- = y^ {sin PLm + cos j3Km)Dm- (26) 


m—O 


where 


D2fc = (0|at2fe|0)„ = (0|at2fc|0)^ = 


1 {2k)\v 


I ,,k 


fu 2^k\ ’ 

D2fc+1 = (l|a^""|l)c = -*(l|at""|l)a = 


_/iut2feh\ -_,7ii7vt2fcn\ _ J_ (2fc-H)! ^ 


3/2 2^k\ ■ 

where k = 0,1,2,.... The results are consistent with the 
isotropic two-photon Rabi model ones, A = lEl. 

We comment that the construction of G-functions as in 
XifM lies ultimately on the fact that the spectral problem 
for the model ([8]) can be recast (through a differential re¬ 
alization of sm(1, 1)) to a differential equation of the Heun 
type. The analytical properties of the G-functions, giving 
in turn the eigenvalues of the model, correspond to spe¬ 
cific conditions of analyticity that Heun functions must 
fulfill to be indeed well behaved solution of the spectral 
differential equationfl^. See [2^ for a neat derivation of 
the G-functions as properties of the Heun functions. 

As for the single-photon Rabi modelfl^, the poles of 
G(z) provide the eigenvalues of the uncoupled bosonic 
mode, A = 0. This is obtained by putting am = 0 in 
Eq.d 




(24) 

(25) 


( 27 ) 
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FIG. 1: (color online) Function of two-photon 

anisotropic Rabi model for lj = 1, A = 0.2, g = 0.3, A = 0.25 
for even (a) odd (b) Fock subspaces, respectively. The blue 
(green) lines are G+ (G!j_), and red (purple) dashed lines are 
G1 (GIL), the poles of have labeled by dotted lines. 


for m = 0,1, respectively. For the parameter in Fig. [2l 
A = 0.2, A = 0.25, we get 50 = ~ 0.6532, = 

~ -0.2483, and ^ ^ 0.3771, = 



I « 0.7689. The positions of level-crossing are 


consistent with them very well. 


IV. CRITICALITY IN THE ENERGY LEVELS 

Despite our model describes a coupling between a 
bosonic field and a two levels system in ‘zero dimension’, 
the energy levels of the system display a ’critical’ behav¬ 
ior resembling very much that one occurring in quantum 
phase transitions of many-body systems. 


A. Discontinuity of entanglement entropy 

First, we discuss the wave function: 


where r]' 


1 - (1 - A) 


2 g 


and m = 0, 2,4,... for 


m = 1,3, 5,... for G^, respectively, as shown in Fig. 
[TJ In the full fledged interacting case, the ‘regular spec¬ 
trum’ are given by the zeros of G^’^ as The ‘ir¬ 

regular spectrum’ providing the well known isolated in- 
tegrability conditions in the parameter space A, g entails 
= 0, for m = 0,1, 2,.... Such points emerge 
as degeneracy point between odd and even sectors. 


4*0 


G^i 


or 


^2 \ 
) 


(29) 


where the constant G = ±1 or G = ±*. In each parity 
sector, 4'c' reads as in (IB19I) . At parity degenerate point, 
the wave function with fixed parity G can be determined 
as follows. 


( Em=Q CLm\m)) + K„,\m)) \ 
V Em=o CKm\m)) + Lm\m)) ) 


(30) 



where M is the truncated number, |m)) denotes \m)a, 
in the appendix it can denote |no -I- m)b in two-mode 
case, such that the wave functions can be written in a 
unified form, similarly |m)) is the corresponding basis 
under transformation > ia^. For our anisotropic 
model, the first energy crossing point corresponds to 
iF^+i(£;r'") =0, m = 0,l: 


Kn = 1 , 


Lq = tan 2/3, 


(31) 


and the corresponding wave function is 

f tan2/3|0)) \ 

' | 0 )) ■ 


(32) 


FIG. 2: (color online) Spectrum of the two-photon anisotropic 
Rabi model for A = 0.2, A = 0.25 in even subspace m = 
0, 2, 4 (a) and odd subspace m = 1, 3 (b). The solid lines are 
the numerical solutions by diagonalizing Hamiltonian, and 
eigenvalues (E^,+) and E’^ _ {E^ _) are labeled by blue 
(green) circle and red (purple) square, respectively. 


Therefore, the wave function of the Z 2 Hamiltonian is 


4 - = 


(cos /3 -h sin/3)|0)) \ 
(cos /3 — sin/3) jo)) J ’ 


(33) 


another eigenfunction is 


The first Juddian solution can be obtained from 60 = 0 
and 61 = 0 , 


A = (1 - A- 


2\ 9 
2w’ 


(l-A^) 


2oj ’ 


(28) 


(cos/3-sin/3)|0)) \ . 3 . 

y (cos/3-b sin/3)|o)) ) 

The eigenfunction labeled by parity G can be written as, 
4-^ = GT'-bl'. (35) 
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g/co 

FIG. 3: Entanglement entropy of the ground and first excited 
state of 2-photon anisotropic Rabi model of even (a) and odd 
(b) Fock subspace of Bogolubov operator, respectively, the 
parameter is the same as Fig. [2l According a level crossing 
occurs at g = (a) and (b) in the energy level, there is 

a sharp discontinuous parity change of the ground state from 
-I- to — as g become larger, as shown by a blue-dot line. 


The entanglement entropy is 

S = —trplogp, (36) 

where p is the reduced density matrix obtained by tracing 
out the bosonic degrees of freedom p = trf,oso„jc(|^')('I'|) 
[ 2 ^. For the first two levels: S = — 
where is the eigenvalue of p. See Fig. O we can find 
that the entanglement entropy of the ground state jumps 
at a critical point. The size of the jump depends on 
anisotropic parameter A. 

We comment that the discontinuity displayed by the 
entanglement occurs with the same mechanism (level 
crossing) in which first order phase transitions occur in 
many body systems psl - l^ . We also comment that the 
parity of the ground state changes at this critical point. 

B. Energy levels condensation 

With a logic that is similar to the one applied to second 
order quantum phase transitions, we analyze the behav¬ 
ior of the energy spectrum of H for spin-boson coupling 
in the neighborhood of \g\/ijJ = Pc = 
of generality, we set 5 > 0, A > 0, and a; = 1). Indeed, 
by inspection of Eg. lfni) . we notice that 77 = 0 at 
This implies that all the levels condense at Eq = — i, 
un-respective of m, therefore with a clear parity symme¬ 
try breaking. The transcendental function in turn, 

displays a pole of higher order. 

This scenario indicates that (a) can be non vanishing at 
g = Pc meaning that the bosonic mode is macroscopically 
occupied. Therefore, such energy condensation implies 
super-radiance. 



n 


FIG. 4; (color online) (a) Energy levels for g > Qc and large 
n. The lowest 12 energy levels for g > gc = 0.8 are displayed 
which seem be grouped into three separate bold curves, each 
bold curve actually includes four lines colored black, green, 
red and blue, from down to up, corresponding to the subspace 
of p = —i, —1, i and 1, respectively. The inset displays the 
results for N = 1000. (b) Photon number probability of the 
ground state for g — 0.85 and N = 200. 


Notwithstanding we cannot analyze the energy for 
g > Pc with our exact solution (because of the topo¬ 
logical constraint coming from the non-compactness of 
su(l, 1) algebra discussed in Sec(|ni)), the numerical data 
shown in Fig. 21 indicate that the low energy levels 
decrease linearly, with a slope depending on the total 
bosonic number as —N^ with x^l.l. If N ^ 00 , 
the curve will be a vertical line. Some insight in the 
limit of large n can be acquired by resorting to the 
RWA.Accordingly, the lowest energy is Eq = (n-|-l)a;/2 — 
a/(w - A)2 -h g‘^[n + l)(n -I- 2) Ri (n-|-l)a;/2-5(n-|-|) Ri 
n{gc — g), in the limit of n ^ 1 (the term of a; — A can be 
neglected in the resonant case). This indicates that for 
g > Pc the ground state energy decreases as Eq r; —nuj 
for large n (Fig. SKa)). 

In the FigHDb) where g = 0.85 which is larger than pc, 
the photon number probability of ground state appears 
to be equally distributed among a large number of modes 
n. Accordingly, the entanglement entropy saturates the 
bound 1, as shown in Fig. [3] 

We remark that such a parity change occurs also in the 
anisotropic Rabi model proposed in Ref. [ 2 ^. 

Similar findings were recently reported for the isotropic 
Rabi model by Plenio and coworkers [s^. 

V. DISCUSSIONS AND SUMMARY 

In this paper we studied a two-boson spin-boson model 
of the Rabi type in which the rotating and counter¬ 
rotating spin-boson couplings act with two different pa¬ 
rameters. The exact solution of the model is pro¬ 
vided through a variation of the Braak method. We 
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remark that our formalism can be straightforwardly ap¬ 
plied to two-modes Rabi models enjoying the same alge¬ 
braic structure sm(1, 1 ) 0 Z 2 of the two-boson case (see 
the appendix). The case with explicit parity symmetry 
breaking (inclusion of a term proportional to cTx) could 
be approached as well within our scheme with minor 
changes}^. 

We comment how the spectrum of the model is found 
modulo certain restriction in the system parameters 
(Eq.JTH)). Ultimately, such a constraint arises because 
the procedure leading to the exact solution involves a di- 
agonalization in sm(1, 1) which is a non-compact Lie al¬ 
gebra. In order to go beyond such restriction, one should 
use the continuous representation of Limblad-Nagel in¬ 
stead of the discrete one(^. To this route, a separate 
study will be dedicated. 

Remarkably, despite the Hamiltonian describes just 
two interacting degrees of freedom, the energy spectrum 
of the system displays characteristic quantum critical fea¬ 
tures. We believe that this is possible because of the 
specihc Z 2 symmetry of the system, acting in a ‘zero di¬ 
mension’ space. Specifically, a ’first order quantum phase 
transition’ occurs because two states belonging to differ¬ 
ent parity sectors with the lowest energies cross. Such a 
level crossing is visible as a discontinuity of the entangle¬ 


ment entropy (obtained tracing out the bosonic degree of 
freedom). The other transition occurs at higher values of 
the spin-boson coupling strength. It is displayed as a par¬ 
ity symmetry breaking transition allowing a macroscopic 
occupation of the bossing mode with a ‘off diagonal or¬ 
der’ (a) ^ 0. By exact analysis of the entanglement en¬ 
tropy we established that the energy condesation occurs 
continuously, sharing similarities with the super-radiance 
phenomenon. Finding the critical properties of the sys¬ 
tem within the exact formalism defines a major direction 
for future investigation. 
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Appendix A: A realization of the anisotropic 
two-photon Rabi model in circuit QED. 

In this section, we provide details on how a two-photon 
Hamiltonian can be realized through circuit QED. 

We refer to the phase qubit, interacting with a mi¬ 
crowave fieldiUlll. For such system, three or more lev¬ 
els participate to the quantum dynamics. Phase qubit 
can be realized with a Josephson junction in a high- 
inductance superconducting loop biased with a flux suf¬ 
ficiently large that the phase across the junction sees a 
potential analogous to that found for the current-biased 
junction. The energy differences between the ground and 
first excited states wqi, and the first and the second ex¬ 
citation uJi 2 are different each other. By tuning the bo¬ 
son field frequency to half of the difference between the 
ground and second excitation energy, the first excitation 
is excluded from the dynamics. The transition between 
the ground and the second excitations can occur by two 


photons. In this regime, the Hamiltonian of the magnetic 
flux-biased phase qubit can be written as 

ffph = cOrO^a + ujqa~^a~ J- g(cos9az — shi9ax){a^ + 

Experiments, nowadays can access the ultra-strong cou¬ 
pling regime. Indeed, typical numbers involved in the 
experiments are woi/Stt = 5.5GHz, uii 2 = 4.5GHz, and 
coupling constant g = O.lGHz. These parameters implies 
a detectable Bloch-Siegert shift ^ 9 ^/A = IMHz. 

Anisotropies in the spin-boson coupling can arise from 
mutual inductance between the phase qubit and the con¬ 
trolling SQUID. In this case the Hamiltonian of the 
circuit is Hph + M{a? — a^^), which, apart from the 
gcos9<7z{o? +a^^), term can be recast to the anisotropic 
two-boson Rabi model Eq. CD- 


Appendix B: The anisotropic 2-mode Rabi model 

The Hamiltonian for the two-mode Rabi model reads 

H 2 -m = a\ai + a\a 2 + z + g{(J^ aia 2 + cr~ a\a\) 

-|-5A((7“''aja2 + (T~aia2). (Bl) 

With two modes in this Hamiltonian, the sm(1, 1) Lie 
algebra is spanned by 

Kq = -{a\ai + a\a 2 + 1), K+ = a\a\, K- = 0102 . 

Here, k = k > 0 for two-mode case, every k 

corresponds an irreducible subspace. These subspace can 
also be written as {[no +n-,n)}, nQ,n = 0,1,2,..., here 
|ni, 712 ) (n-i > 712 ) is the Fock state of the boson operator 
ai, 02 . We note that ni < 712 give a set of equivalent 
subspace, {|n, no+7i)}, giving degenerate eigenvalues and 
eigenfunctions; in the following we will work with the first 
choice. Thus we can also recognize tiq as the subspace 
index, the relation of uq and k is uq = 2k — 1. 

The Hamiltonian enjoys the Z 2 symmetry generated 
by 

n2_^ = ^ cos(7ra^a2) (B2) 

which have value p = ±1 in the corresponding subspaces 
of K. 


1 . Exact solution of the anisotropic 2-mode Rabi 
model 

The exact analysis of the two mode Rabi model (IBII) 
proceed along very similar lines we followed for the two- 
boson case. We apply U and the Bogolubov transforma¬ 
tion: 
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- u;ala2 + p H 
—q + sa\a\ 


r(aia2 
f taia2 


a\a\) —q + saia 2 

uja\ai + uja\a2 — p — r{aia2 + a\al) 


ta\a\ 


(B3) 


bi = uai + va^, 62 = ua 2 + va\, (B4) 


with 


we can get 


l-7l< 


2uj 


1 + A|’ 


(B5) 


u = 


'1 + C 


2 C 


v = 


cos 2 / 3 =l^C, (B6) 


where 


c = 


'l-(l + A)252/4a;2 


(B7) 


1 — (1 — j^uP- ’ 

The vacuum state can be found in 6 i 62 |uoj 0)& = 0, as 

^ 00 

n=0 


- 00 


+ n, n)(B8) 

V no\n\ 


which also have property 6 i|no, 0 )b = ^/^\no — l, 0 )h. 
Here, is the normalized coefficient, which can be 

found with the Taylor expansion. 


1 

(1 — 


1 + (no + l)a; + 


(no + 2) (no + 1) 
- ^ + 


where x = v^ju^. Actually, |no,0);, is a two-mode 
squeezed vacuum state, which can be generated by the 
2 -mode squeeze operator S' 2 (C) 01 ^ the vacuum 


|no, 0)b = S2{0\n-o, 0) = exp (^^44 “ Caia2^ |no, 0), 

where u = and v = —In quantum optics 

S' 2 (^) is associated with degenerate parametric amplifi¬ 
cation [I|. 

The Fock states can be constructed as 


|no-k TO,m)& = (6j6^)™|no,0)b, (B9) 


Note that, the normalized coefficient y elim¬ 

inated for the simplicity. 

The ansatz for the eigenstates of {H) is 


</>i 

4>2 


Em=o I^m\no + m, ■m)b \ 
Em=o I^m\no + m, m)b ) 


(BIO) 


Because of the Z 2 symmetry, the constant C has only 
two values ±1. So we can construct the transcendental 
function as: 


(Bll) 


use the coefficients above we can obtain G-function as: 

00 

(- cosl3Lm + s\nl3K^)Dm, 

m—0 
00 

<^4- = Y,{sinPLm +cos PKm)Drn. (B12) 

m=0 


where the coefficients Dm is as follows 


Dm = (no,0|(6l6^)™|no,0)6 = (no,0|(44)™l”o,0)^ 
1 /UN ™ (riQ -I- m)\ 


u'^ 0+1 Vm 


© 


no! 


(B13) 



FIG. 5: (color online) Function of two-mode anisotropic 
Rabi model for a; = 1, A = 0.2, g = 0.5, A = 0.5 for no = 0 
(a) no = 1 (b), respectively. The blue lines are G"“, and red 
dashed lines are , the poles of have labeled by dotted 
lines. 


The analytical property of is similar to it in the 
isotropic case, as shown in Fig. [5] The poles of can 
be found similar as that for two-photon case, 

El°m = + ‘2rn + l) - uj. (B14) 


where C 


1 - (1 - A) 


2 5 
4a; ^ 


The energy spectra has 


divided into two parts, one is the regular case, which 
correspond to the zeros of G"“j_ as ^. The other is 
the irregular case, which correspond to the level-crossing 
points, i.e., Juddian solutions. They can be found sim¬ 
ilarly as Km+i{E^°m) ~ then the numerator and de¬ 
nominator of are both 0. However, there is a few 
difference, in anisotropic case the first degenerate point 
is at m = 0 , not m = 1 in isotropic case, because of the 
parameter A of the counter-rotating term. 


^x°± ~ (’2'0i 0|[‘/?2 T 
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FIG. 6: (color online) E'^°^ vs. g of no = 0, m = 0,1, 2 
(a) and no = 1, m = 0, 1 (b) for a; = 1, A = 0.2, A = 0.5. 
The solid lines are the numerical solution, and eigenvalues 
of parity + and — are labeled by blue circle and red square, 
respectively. 


The Juddian solution at m = 0 can be found as 


So the eigenfunction can be written as 



or 


( ^o|0)) \ 

I L,m) 


(B18) 


where basis |0)) can be |no:0)6: |0)a and |l)o,, and |0)) is 
corresponding case of |0)). 

Using parity they can be written as 




( CLo\0))+Ko\0))\ 
\CKo\0))+Lo\0)) j ’ 


(B19) 


where C = ±1, ±f, the two eigenfunctions are orthogonal 
to each other. 


We can study the entanglement entropy of the ground 
state. Similar as in two-photon case, here we can also 
find that the entanglement entropy of the ground state 
jumps at a critical point, see Fig. 0. 


+ (B15) 

if A is known, from bo = 0 the especial value of g can be 
found as 

As shown in Fig. 1 we can exactly determine the 
first level-crossing points at subspace ng = 0,1 for pa¬ 
rameter A = 0.2, A = 0.5. In no = 0 subspace, it is 

g = ^ ^ 1.0328, E = yil - 1 « -0.3890. And 

g = ^ Ki 0.7303, E = - 1 Ri 0.6452 for no = 1 

case. The wavefunctions at these degenerate points will 
be discussed in detail later. 


2. Entanglement entropy 

For example, in isotropic 2-mode and 2-photon case, 
the first level crossing occurs in m = 1, 

A 

2w?7 



FIG. 7: Entanglement entropy of the ground and first excited 
state of 2-mode anisotropic Rabi model of no = 0 (a) and 
no = 1 (b), respectively, the parameter is the same as Fig. [6] 
The level crossings occur at g = (a) and (b) in the 

energy level, there is a sharp discontinuous parity change of 
the ground state from + to — as (? become larger, as shown 
by a blue-dot line. 


Ko = 1 


(B17) 





























